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We study a composite milliciiarged dark matter model. The dark matter is in the form of pion- 
like objects emerging from a higher scale QCD-like theory. We present two distinct possibilities 
with interesting phenomenological consequences based on the choice of the parameters. In the first 
one, the dark matter is produced non-thermally and it could potentially account for the 130 GeV 
Fermi photon line via decays of the "dark pions". We estimate the self- interaction cross section 
which might play an important role both in changing the dark matter halo profile at the center 
of the galaxy and in making the dark matter warmer. In the second version the dark matter is 
produced via the freeze-in mechanism. Finally we impose all possible astrophysical, cosmological 
and experimental constraints. We study in detail generic constraints on millicharged dark matter 
that can arise from anomalous isotope searches of different elements and we show why constraints 
based on direct searches from underground detectors are not generally valid. 

PACS numbers: 95.35.-|-d 95.30.Cq 



I. INTRODUCTION 

The current knowledge is that all particles in the Stan- 
dard Model are either electrically neutral (e.g. neutrinos, 
photons and gluons), or have charges which are integer 
multiples of e/3 where e is the charge of the electron. The 
possibility of particles with electric charges much smaller 
than e/3 has always intrigued both theorists and exper- 
imentalists. One way one can add millicharged particles 
in the Standard Model is for example by adding a parti- 
cle with hypercharge F = 2e that is singlet of the color 
SU(Z) and weak SU{2). However as it was pointed out 
in [1[ , the embedment of such an extension of the Stan- 
dard Model into a Grand Unified Theory is problematic. 
These problems can be avoided if an extra f/(l)' is intro- 
duced [2] . Due to mixing between photons and the para- 
photons of the [/(I)', particles charged under the U{1)' 
appear to have a small coupling to the photon and thus as 
electrically millicharged. A third way millicharged par- 
ticles might appear is in the form of Standard Model 
neutrinos that have tiny electric charges. A small modi- 
fication of the hypercharge assignments of the Standard 
Model particles can accommodate tiny electric charges 
to neutrinos without inducing any gauge anomaly (for a 
review see e.g. 0-Q)- 

Although Charged Massive Particles (ChaMPs) have 
been more or less ruled out as dark matter candidates, 
this is not easily the case for particles with tiny charges. 
Millicharged particles have been firstly proposed as a 
solution of the dark matter puzzle long time ago [6|- 
Q. Composite dark matter (with charged constituents) 
is an attractive scenario that has been studied exten- 
sively [9I-KL4J. Recently there has been a revived interest 
on millicharged dark matter via Stueckelberg Z' mod- 
els flSl Hq , in the form of millicharged atomic dark mat- 
ter [l7[, or in the form of MeV particles in order to ex- 



plain the 511 keV line [18|. 

In this paper we present a model of composite dark 
matter inspired by QCD. As in QCD where the sponta- 
neous symmetry breaking of the chiral symmetry leads to 
Goldstone bosons, similarly here, chiral symmetry break- 
ing produces Goldstone bosons that play the role of dark 
matter. We explore the phenomenology of this theory ex- 
amining two interesting cases. In the first one, the "dark 
pions" are produced non-thermally and they can explain 
the 130 GeV Fermi photon line via decays of the form 
ttq — > 27. In the second version, we show how the pions 
can be produced thermally via the freeze-in mechanism. 
We also study all possible astrophysical, experimental 
and cosmological constraints imposed on such a model. 

We chose a QCD-type Lagrangian of the form 



L — Tpi^f^D'^^p — mip^p ■ 
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where the gauge group is an SU{3) "color". We as- 
sume two flavors of fermions that carry tiny electric 
charges 2ee/3 and — ee/3 where e is a small number that 
will be estimated later on. The covariant derivative is 
Dt' = d'' - igr'^G'^'' - iqA'', where G"^ is the "new 
gluon" field and A^^ is the electric field. We assume that 
the theory becomes confined below a scale, where the de- 
scription of the theory in terms of mesonic and baryonic 
degrees of freedom is more appropriate. There are three 
Goldstone bosons i.e. the pions of the theory 7r+, 7r~, and 
ttq. As in QCD, it is expected (due to electromagnetic 
interactions) that n are heavier than ttq. For simplicity 
we are going to assume that there are no baryons of the 
theory present, and therefore dark matter is composed of 
the three pions of the theory. This is expected in cases 
where the "new baryon" number is not conserved. We 
study now the phenomenological implications of this sim- 
ple QCD inspired dark matter model, starting from the 
possibility of producing the Fermi 130 GeV line. 



II. THE 130 GEV LINE 

Recently, analysis of the Fermi collaboration data re- 
vealed the possibility of the existence of a Tline at 130 
GeV coming from the center of the galaxy |19l-l24l| . Due 
to the fact that this line is sufficiently peaked and does 
not follow the typical astrophysical background power 
law behavior, this can be an indirect signature of dark 
matter presence. In general, if dark matter is in the form 
of WIMPs (denoted by x), processes like X + X — >■ 7 + 7, 
X + X~^ l + Xj or X ^ 7 + 7 are in principle all allowed. 
Obviously, in the case of annihilating dark matter to two 
photons, the WIMP mass must be 130 GeV, and in the 
decaying dark matter scenario, the mass must be 260 
GeV. For the case where one photon and one other par- 
ticle X is produced, the WIMP mass depends obviously 
on the mass of X. One striking observation is that if the 
line is attributed to annihilation of WIMPs, the annihi- 
lation cross section does not seem to be the one required 
for thermal production of WIMPs, but it is rather smaller 
by roughly an order of magnitude depending on the as- 
sumed dark matter halo profile. This could be an indi- 
cation towards an asymmetric WIMP annihilation where 
dark matter and anti-dark matter do not appear in na- 
ture with equal numbers (something supported by some 
recent analysis |25[) or towards a decaying dark matter 
scenario. However, in both cases a large branc hing ra - 
tio to photons (larger than ^ 0.01) is required J26l - l28l |. 
Otherwise, decays and/or annihilations to other parti- 
cles can induce a continuum spectrum of photons that 
can hide the line. From this perspective a dark pion that 
has only the decay channel ttq — ^ 27 fulfills perfectly this 
requirement. 

Although Einasto or Navarro-Frenk- White dark mat- 
ter halo profiles favor the scenario of annihilating dark 
matter [19j, decaying dark matter could account for the 
130 GeV line if the dark matter morphology at the cen- 
ter of the galaxy is appropriate (i.e. more cuspy). The 
photon spectrum of a decaying WIMP to two photons is 



d$H 



dE^dfl in 



re 



Pe 
2m^ 



1 f Phe.\o{r)\ dNd. 



ds . 

los '"O V PO 



dE^ 



while for annihilating one 



(2) 



rf$a 



dE.ydn 



{av} I PQ_\ I , J_ f Phe,io{r)\ dN, 



Stt 



my 



los 



^■0 V p& 



dE~, 



where P, and {av) are the decay rate and annihila- 
tion cross section to two photons respectively, dN/dE^ 
denotes differential photon energy spectrum, pq ~ 
0.3 GeV/cm'^ is the local dark matter density, Tq ~ 
8.5 kpc is our distance from the center of the galaxy and 
Phaio('^) is the dark matter halo profile. The integral is 
taken along the line of sight. Decaying dark matter has 
already been proposed as an explanation of the 130 GeV 
line m^- 



In this paper we do not attempt to make a sophisti- 
cated analysis of the Fermi spectrum. Following [2^ on 
generic grounds, we assume that the observed spectrum 
is attributed to decaying dark matter if one adjusts the 
decaying rate such that the spectrum will be similar to 
the annihilating one. Depending on the chosen halo pro- 



file, an approximate value of P 
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satisfy the requirement [29| . We should mention that the 
analysis in [1^ suggests that annihilating dark matter 
fits much better the Fermi spectrum compared to a de- 
caying dark matter scenario. However a decaying dark 
matter scenario can be valid if the dark matter density is 
more cuspy at the galactic center. Although this is not 
currently favored by numerical simulations of collision- 
less dark matter, the existence of self-interactions might 
alter completely the situation. We shall show that indeed 
such self-interactions are present among the dark pions. 
The ttq of the theory decays due to the anomaly exactly 
as in QCD. However, due to the fact that the constituent 
fermions have tiny charges, the lifetime time of the par- 
ticle can change significantly. The decay rate is 



a'^e^ ml^ 



647r3 /2 



(4) 



where a is the usual fine structure constant. Demanding 
the decay rate to be ~ 10"^^ sec~^ gives the condition 



6 X 10-l^ 



(5) 



where /^ is measured in GeV. For a typical value of 
/^ = 1 TeV, e ~ 1.9 X 10"^^. The pions can also co- 
annihilate. In order to find what is the annihilation cross 
section, we should introduce an effective chiral perturba- 
tion Lagrangian that describes the mesonic fields. Simi- 
larly to QCD 



f2 j-2 

L = ^Ty{D^J:D''T.'') + ^Tr( 
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where D^ — dfj_ + ieA^[Q, E], Q being a diagonal matrix 
with entries the electric charges of the two "quarks" in 
units of e. The matrix S is defined as usual 
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where r" are the Pauli matrices. It is easy to see that to 
lowest order ttq does not couple to photons (as it is ex- 
pected), because the commutator [Q, E] is zero. However, 
ttq's can annihilate to photons through one-loop media- 
tion of TT^ states that couple to photons. There are two 
Feynman diagrams contributing to the amplitude. For 
on-shell production of photons the total amplitude is 

M = M^.e'^e^-', (8) 

where e^'^ represent the helicity of the two produced pho- 



tons and M^v is [32 
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where zj- = 1 ± (1 — 4?nJ/s)^'2, k represents photon mo- 
mentum, and s is the typical Mandelstam variable. It 
is understood that for the derivation of this result, we 
have made the approximation that the mass difference 
between tt and ttq is not large. Given the above ampli- 
tude for ttqtto — >■ 77, we can compute the corresponding 
cross section in the center of mass frame 
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where /(s) = 2[\n^{z+/z^) - tt^] + {ml/s)[\n'^{z+/z^) + 
Tr^]^. Upon making the assumption that ttq move with 
non-relativistic velocities (s ~ 4m^ and z± ~ 1), and by 
using Eq. (0, 
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This annihilation cross section which in principle can pro- 
duce a photon line at 260 GeV is so tiny to be detected 
or to play any role in the evolution of the dark matter 
density. 

It is easy to show that tt's are not in thermal equilib- 
rium with the plasma for the e values we are interested. 
We can have a rough estimate of the Thomson cross sec- 
tion between a charged pion and a photon. The condition 
for thermal equilibrium is that n^axv > H , where n-y is 
the number density of photons. The Thomson cross sec- 
tion (in natural units) is ctt = {8TT/3in'^)a^e'^. Using 
rij = (2C(3)/7r2)T3 and H ^ IMgl^'^T^/nipu it is easy 
to show that the dark pions (or dark fermions before 
confinement) with e < 10~^ were never in thermal equi- 
librium with the plasma (unless one goes to extremely 
high temperature when pions did not exist). This means 
that in this case tt's have to be produced non-thermally 
by some unknown mechanism. Whatever the mechanism 
is, it should not disturb the thermal evolution of the uni- 
verse (i.e. the effective degrees of freedom g* should not 
change). This can easily happen for example if the pi- 
ons are produced via some decay with relativistic ener- 
gies but with a temperature significantly lower than that 
of the plasma. As we pointed out above, the plasma is 
not in thermal equilibrium with the pions and therefore 
they do not have to have the same temperature. How- 
ever pion self-interactions are sufHcient to establish ther- 
mal equilibrium among the pions providing pions with 
a common temperature that must be much smaller than 
the temperature of the plasma in order to avoid prob- 
lems with the evolution of the universe. The reason that 
self-interactions can keep pions in thermal equilibrium is 



due to the fact that the cross section is independent of 
e and therefore is not suppressed. The coupling among 
pions can be estimated from Eq. ^ if one expands the 
second term to fourth order in the fields. The amplitude 
for ttqttq •!-> TToTTo (but also for tt+tt^ — > ttotto if one as- 
sumes that the difference in the masses of ttq and tt^ is 
small) is 
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This leads to a cross section 
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The condition riT^aTr-nV » H is easily fulfilled. One 
can check by using the number density for a relativis- 
tic species with a temperature T^m << T~^ and the above 
cross section, that pions were in thermal equilibrium with 
each other at some point. Obviously as the universe ex- 
pands and the pions become non-relativistic, pions at 
some point will lose contact with each other, but as in 
the case with neutrinos, they will preserve their Bose- 
Einstein distribution since their entropy is preserved. 
The fact that pions were in thermal equilibrium enable 
us to estimate how much tt^ is left today. Assuming that 
TT* are slightly heavier than ttq (as it happens in QCD), 
the annihilation cross section for 7r"'"7r~ — >■ ttottq is given 
again by Eq. (J13p . The difference between the annihi- 
lation and the self-interactions among pions is that the 

square root in Eq. (jl3|) is 



'1 — rri^^Jm^ for the for- 
mer and ^ V for the latter. It is understood that in the 
case of annihilation at the nonrelativistc limit s ~ Am? + . 
In the non-relativistic limit, the annihilation cross section 
takes the form 



2.3 X 10-^^cm' 
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where /i is the value of f-„ in units of 1 TeV, and m27o 
is m^± in units of 270 GeV. This is a quite interesting 
result. By playing with /^ and the mass splitting between 
TT^ and ttq, one can get from complete extinction of tt* 
today to a partial one. In the case of partial presence 
of TV , assuming that x fraction of dark matter is in the 
form of ttq and 1 — a; in the form of tt , Eq. ([5]) should 
be modified to 
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This is because the Fermi 130 GeV line has to come from 
decay of the x part of the dark matter density (i.e. ttq). 
There are two comments in order here. Due to the small- 
ness of e, no essential constraints can be imposed on the 
1 — a; part of dark matter that is made of tt* (see the 
constraint section). 



The second comment has to do with the fact that the 
freeze-out of the annihilation takes place at a tempera- 
ture Tdrn which as we shall show below must be signif- 
icantly smaller than the plasma temperature T^. How- 
ever, since everything takes place in high enough energies, 
g^, and f;,s do not change substantially, and therefore one 
can approximately consider that the freeze-out occurs as 
usually at T^ai ~ m^±/25. However one should keep in 
mind that this corresponds to a much higher plasma tem- 
perature. One would naively expect that earlier freeze- 
out would correspond to even colder dark matter. This is 
not necessarily true due to self-interactions. We discuss 
this issue later in this section. 

As we mentioned Tdm is in principle different than T^. 
If pions have been produced relativistically, one has to 
make sure that by the time of the Big Bang Nucleosyn- 
thesis (BBN), dark matter does not change significantly 
the effective number of relativistic degrees of freedom g^ . 
In order to find how the constraint affects our dark mat- 
ter scenario, let us assume that the pions are produced by 
a time when their temperature was T^i and the plasma 
temperature was Ti. Let us call /? = Tdi/Ti. Now let's 
assume that the dark pions become non-relativistic at 
some time where their temperature is Td2 — in and the 
plasma's temperature is T2. Up to that point both the 
plasma and the pions are relativistic and therefore the 
ratio T(i2/T2 is still /3. Below that temperature pions be- 
come non-relativistic. However, since pions were never 
in thermal equilibrium with the plasma, their entropy 
should be conserved. At the BBN time scale, entropy 
conservation dictates riBBNOBBN = '^2«2 (where 2 refers 
to the time where pions become non-relativistic). The 
energy density of pions at the BBN time scale is pbbn = 
TO^BBN = ?^2™(a2/aBBN)^ = ?^2™(T'BBN/^2)^, where 
Tbbn = 1 MeV. Using the fact that na = C,{'i)Tl^/-K'^ 
and that T^^jT^ = /?, one can derive the contribution of 
pions to (7* at BBN as 



the dark matter density today. This leads to 
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Extra relativistic degrees of freedom at the epoch of BBN 
are tightly constrained because the abundances of light 
elements are strongly affected by 5*. At BBN g^ = 5.5 -I- 
7/4 * Ni^ where Ni, is the number of neutrino families. 
Analysis of the BBN light element abundances and CMB 
data from Planck limit N^, down to ~ 0.5 [33| (on top 
of the three existing neutrino families). This leads to a 
freedom of gX' < (7/4)0.5. This dictates that /3 < 0.02 
in order not to destroy the BBN predictions. We can now 
estimate what is the /3 needed in order to get the observed 
relic density of dark matter today. In order to do this, 
we use again the conservation of the dark matter entropy 
"bbnObbn — '^O: where ng = 1.176 x lO^^GeV/cm"^ is 
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The only unknown cbbn can be easily estimated by using 
conservation of the entropy for the plasma this time i.e. 
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gl.T^. Using g 



BBN 



10.75, gl. 



3.909 and Tq = 2.728 K, we can get obbn = 1-68 x 10"^°. 
Using Eq. (IT7)) leads to a value /3 ~ 4 x 10^** which is 
much smaller than the 0.02 constraint we derived above. 
As we have already mentioned, self-interactions are in- 
dependent of e and are not suppressed by powers of it. 
This can be quite crucial for two reasons: Firstly self- 
interactions might change the dark matter halo profile at 
the center of the galaxy. Although with standard dark 
matter halo profiles, Fermi data (19{ favor annihilating 
dark matter over decaying, dark matter self-interactions 
might change completely this picture, since most dark 
matter simulations so far consider collisionless dark mat- 
ter. The second point is that self-interactions might 
make the dark pion to be a warmer dark matter can- 
didate which is supported by some data [3J|. The non- 
relativistic self-interaction cross section can be read off 
from Eq. ^ 
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The value is sufficiently small to enable the dark pion 
evade constraints on dark matter self-interactions im- 
posed by the bullet cluster [35L, or galaxy elliptic- 
ity [3a, I33], or old neutron stars 38l-l40|. However, the 
cross section being of the order of picopbarn or even 
higher (if one selects a smaller /^r and/or larger TOtto), can 
potentially have an effect on how warm the dark matter 
candidate is, and on the inner dark matter halo profile. 
The answer to especially the second issue is difficult and 
in principle requires simulations of self-interacting dark 
matter which might be considerable different from the 
collisionless case. 



III. THERMALLY PRODUCED 
MILLICHARGED DARK MATTER 

Despite the tiny electric charges of its constituents, 
and the fact that tt's are not in thermal equilibrium with 
the plasma for small values of e, they can be produced 
thermally via the freeze-in mechanism [41| . The process 
77 — >■ TT^Tf- is sufficient to produce an abundance similar 
to the dark matter one. The Boltzmann equation for the 
process above can be written as [4l| 
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It is not hard to show that the leading s contribution 
of the cross section conies from the Feynman diagram 
with the quartic couphng e^e^A^A'^7r+7r~. This diagram 
contributes to |AfP _^^+^_ = 2e'*e'*. Integrating Eq. p^ 
gives 
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where x — vrij^±lT^ and s\ — (27r^/45)(7*sT'^ is the total 

entropy. Using t ~ O.3OI5* (mpi/rn^±)a:^, the equa- 
tion can be rewritten as 
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We can now integrate the equation from a; = to infinity. 
Using 
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and g* ~ g^a ~ 100, we get e — 1.3 x 10^^. This means 
that with e ~ 10~^, the abundance of tt^ is that of the 
dark matter. From Eq. ([5]), one can realize that mil- 
licharged dark matter could in principle be produced via 
the freeze-in mechanism and provide via decays a ~ 130 
GeV line for /^ ~ 4.8 x 10^^ GeV which is the GUT scale. 
However, this is not possible. By a simple rescaling of 
Eq. ([13]), one can see that the value of /^r is so large (cor- 
responding to a tiny tt+tt^ — > ttottq cross section) that 
makes impossible to create substantial quantities of ttq. 
In addition, we found that the direct production of ttq 
from the process 77 — )■ ttotto is suppressed parametrically 
by ("^TTo/Zir)^- This can be seen by using Eq. © and fol- 
low the steps leading to Eq. (PT|) keeping only the leading 
term in s. This means again that no ttq are produced and 
therefore the Fermi 130 GeV line cannot be explained by 
this version of millicharged dark matter. However, as we 
demonstrate in the next section, this thermally produced 
(via the freeze-in mechanism) millicharged dark matter 
evades all possible constraints arising from dark matter 
direct searches in underground detectors. 



IV. CONSTRAINTS ON MILLICHARGED 
DARK MATTER 

The presence even of tiny electric charges on particles 
(and in particular dark matter ones) in nature can lead 



to severe constraints. This has been addressed long time 
ago in the context of ChaMPs but it was extended in 
the case of fractionally charged particles [a-|8|, 14^] . These 
bounds have been obtained by a variety of different meth- 
ods ranging from accelerator searches and dark matter 
searches to constraints from the Lamb shift. Big Bang 
Nucleosynthesis, the Supernova 1987A, plasmon decay 
in red giants and white dwarfs. In addition there are 
cosmological bounds based on the fact that $7 < 1. A 
summary of these bounds can be seen for example in |5| . 
All these bounds are avoided by both versions of the dark 
pion we analyze in this paper. In the 260 GeV decaying 
pion, e is too small to be constrained by these bounds. 
The same is true for the thermally produced dark pions 
which although have e ~ 10~^, it is still too small to be 
excluded (especially if pions are embedded in a model 
with paraphotons). 

In addition, there are severe constraints on e based on 
the bullet cluster [3a| and the ellipticity of galaxies like 
NGC 720 [21 113 ■ The bullet cluster imposes an upper 
bound on dark matter self-interactions and since charged 
particles with charge ee can exchange photons, an upper 
bound is set for e. Coulomb-type self-interactions can 
also make the distribution of dark matter more isotropic 
and this leads to more spherical galaxies. The ellipticity 
of certain galaxies like NGC 720 imposes again an upper 
bound on e. Furthermore, millicharged particles can in 
fact cause changes in the Cosmic Microwave Background 
spectrum leading to very strict constraints on e |43l . |4J| . 
Among the above three constraints, the CMB ones are 
the most strict. However, a millicharged particle with 
mass larger than ^ 100 GeV and e < 10~^ clearly avoids 
it. Below we will study in more detail two other sets 
of constraints i.e. the ones from direct search detection 
in underground detectors and the ones from searches of 
anomalous isotopes of different elements. 



A. Direct Detection 

Direct detection search experiments like Xenon [45 
and CDMS ^] have imposed severe constraints on the 
cross section between dark matter particles and nuclei of 
detectors that have been planted in deep sites in order to 
be shielded from unnecessary background radiation and 
thus false signals. Millicharged dark matter can exchange 
photons with the nuclei and therefore constraints on e can 
arise. Such constraints have been studied in the context 
of mirror dark matter [47| and in general [48] . In principle 
such a constraint might be quite severe. Millicharged par- 
ticles scatter off charged nuclei with a Rutherford cross 



section. Of course sufficiently small values of e lead to 
little scattering and thus particles will avoid detection. 
Based on the CDMS data, the analysis of \48\ provided 
an upper bound e ~ 10~* for millicharged dark matter 
with mass above 10 GeV. This limit is of course evaded 
by our decaying dark pion with an e '^ 10~^^, but it 
could exclude the thermally produced dark pion that re- 
quires e ~ 10~^. However, we show that the excluded 
e-region of [4^ is not valid. We shall demonstrate that 
for intermediate values of e (small but not too small), 
millicharged particles lose most of their kinetic energy 
through electromagnetic interactions with atoms before 
reaching the detectors while traveling underground, lead- 
ing to Rutherford scattering with recoil energies way be- 
low the threshold of the detectors thus invalidating this 
type of constraints. 

A non-relativistic millicharged particle can interact 
and lose energy via collisions with electrons according 
to the Bohr-Bcthe-Bloch formula 
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where n — N aZ p j M^ is the electron density of the tar- 
get (Na, Z, and M„ being respectively the Avogadro 
number, the atomic number of the material, and the mo- 
lar mass constant), nif^ is the mass of the electron, and 
/ ~ lOZeV is the mean ionization potential. The energy 
loss above can be trivially rewritten as 



'dx ^ E '^Eb' 



(25) 



where B — 2Trne^e*{m/me) {m being the mass of the 
millicharged particle) and Eq = ml /{Ame). The equa- 
tion can be now integrated from an initial value of kinetic 
energy E{-c^ (corresponding to the usual dark matter ve- 
locity of 230 km/s) to a final Ef. We can easily find 
the upper value of Ef in order for a collision to produce 
enough recoil to trigger the detector. The recoil energy 
T takes values between < T < AEfmniM /{ra + mpf)'^. 
A uniform distribution of T is usually assumed within 
this range. If we demand the recoil energy to be less 
than 1 keV (which is already below the thresholds of 
current experiments), we can estimate what minimum 
value of Ef can produce this recoil. By solving the dif- 
ferential Eq. (|25p and tracking the solution from Ein to 
Ef we can find what distance is needed to slow down 
the millicharged particle to an Ef that cannot produce a 
detectable recoil in the underground detectors. We find 
that for an e = 10~^ and m = 1 TeV, this distance is be- 
low 2 meters. For this we use an average matter density 
of S.Sg/cm'^. The dependence of this result on Z is ex- 
tremely weak and therefore not relevant. CDMS lies 2341 
feet below the surface and Xenon lies below 3700 meters 
water equivalent. Clearly no constraints can be imposed 
on both versions of dark pion: the first one because the 
particle is neutral with constituents of e '-^ 10~^^ that 
is too small to produce a significant rate, and the sec- 
ond (with e ~ 10~^) because the particle slows down 



significantly before it reaches the detector. Assuming a 
depth of 500 meters (which is a bit smaller that both the 
depths of the sites of Xenon and CDMS), we find that 
a TeV millicharged particle can slow down sufficiently 
to avoid direct detection if e > 5.4 x 10^^. Therefore 
the overall excluded parameter space for that mass is 
10~® ^ e < 5.4 X lO^'^. For 100 GeV milhcharged parti- 
cle this range reduces to 10~* ^ e ^ 1-7 x 10~^. We have 
used an "effective" Z = 40 for the ground material. The 
dependence of the above limits on Z is weak. 

Finally there is another point that enables a thermally 
produced dark pion to evade detection. As it was pointed 
out in 49], for lOOe^ ^ 'ti < lO^e TeV, galactic magnetic 
fields prevent millicharged dark matter from the halo to 
enter the galactic disk and those initially trapped in the 
disc are accelerated and injected out of the disc within 
0.1 to 1 billion years. Thus no constraints on millicharged 
dark matter can be imposed based on the direct search 
experiment within the above range. For e ~ lO^"" the 
GeV- TeV mass range of the millicharged particle is inside 
the above non-detection region. 



B. Anomalous Isotope Constraints 

In this section we present generic constraints valid for 
any type of millicharged dark matter particles based on 
searches of anomalous isotopes of hydrogen and oxygen 
in the sea water and helium in the atmosphere. 

The binding energy between a (negatively) mil- 
licharged heavy WIMP with a nucleus is 



£^Rin — 



e^Z^a^mAT 



(26) 



where we assumed that mj\j « m. Typical examples 
of couplings with hydrogen and oxygen nuclei give corre- 
spondingly 



Eh ^ 2.5 X lO^^eV 
Eo c^ 2.6 X lO^e^eV. 



(27) 



The kinetic energy of a WIMP assuming a velocity v 
230km/sec is 



El. ~ 2.9 X 10^ 



GeV 



eV. 



(28) 



First we can calculate the cross section for the cap- 
ture of a millicharged particle by a nucleus. If the Bohr 
radius of a bound state between the nucleus and the milli- 
charged particle is smaller than the corresponding atomic 
one (nucleus-electron system), one can use to good ap- 
proximation the formula that gives the cross section for 
a capture of an electron by a charged ion [50| 



aN 



75„3„2 

-Z a flg 



E" 

E 



5/2 



(29) 



This is the cross section for the capture of an electron 
by a nucleus of charge Ze, where a is the fine structure 



constant, E is the energy of the electron, and uq and 
^Bin ^^^ respectively the Bohr radius and the binding 
energy of the hydrogen atom. This formula has been 
derived at the limit where a moving light electron gets 
trapped by a heavy nucleus at rest. In order to be able 
to use this formula in our case, we have to go to the rest 
frame of the heavy millicharged particle. Therefore in 
our case of interest (i.e. the capture of a milli-charged 
particle with charge — ee by a nucleus with charge Ze), 
Z has to be substituted by Ze in Eq. (f^ . gq is the Bohr 
radius having substituted the mass of the electron with 
the mass of the nucleus since the reduced mass of two 
unequal masses is given approximately by the smaller 
one which in this case is the mass of the nucleus, and 
i?gi,j is the binding energy of the system nucleus- WIMP 
having normalized Z = e = 1. Due to the strong depen- 
dence of capture cross section on the energy E, there are 
two distinct cases here. In the first one, the millicharged 
particle loses enough energy via electromagnetic interac- 
tions with electrons (Eq. ([M])) and slows down to ther- 
mal velocities. Since the millicharged particle is heav- 
ier than the nuclei of the water and air that we are go- 
ing to consider, the particle might be considered at rest 
and therefore E in Eq. 1^^ is the thermal energy of the 
atoms 3kBT/2. If electromagnetic interactions do not 
slow down sufficiently the millicharged particle, E is the 
energy of an atom moving with the velocity of dark mat- 
ter V ~ 230km/sec. 

There are two main constraints from anomalous iso- 
tope searches that we consider, i.e. the hydrogen and 
oxygen ones from ocean water and the helium ones from 
the atmosphere. Using Eq. ((24)) in the case of water, 
we find that within the average 2.6 km depth of the 
oceans, millicharged particles of 100 GeV (1 TeV) mass 
with e > 2.1 X lO"'^ (e > 6.7 x lO"'^) slow down to their 
thermal velocities. Using this fact and Eqs. (|27p and p9|. 
we find the capture cross sections for this regime 



o-H 



1.4 X lO^^e^cm^ 



ao ^ 1.9 X lO^^e^cm^. 



(30) 



The mean free path of a millicharged particle before 
it gets captured by a hydrogen (oxygen) nucleus will 
be given by Xh^q) = l/"-ff(O)0-ff(O), where uh = 
Na/{9 cm^) and no — Na/{18 cm^) are the number 
densities of hydrogen and oxygen respectively in water. 
The probability for capturing a millicharged particle is 
Ph(o) = 2.6 km/A//(o) if Xh(o) > 2.6 km, or one oth- 
erwise. As we mentioned, 2.6 km is the average ocean 
depth. It is easy to check that due to the fact that 
ao >> <Jh , most captured millicharged particles are cap- 
tured by the oxygen nucleus. The ratio of the number 
density of the millicharged particles over the density of 
hydrogen in the sea water will be 



"-X 



Px 



pv 



■U 



n[{m^ 2.6km 



(31) 



where p^ = 0.3 GeV/cm'^ is the local dark matter den- 
sity in the earth, p is the probability of capture derived 



above, v — 230 km/sec is the velocity of the millicharged 
particle, and tacc = 3 x lO^years is the approximate age 
of the oceans. If e is sufficiently large, then the reactions 



H + e^ He 
O + e^Oe, 



(32) 



are in thermal equilibrium. We take the condition for 
this to be itherm = 1/nav << iacc: where v = y^3T/m 
is the thermal velocity of a hydrogen or oxygen atom. 
The condition for thermalization is fulfilled for e > 5 x 
10~^ in the case of oxygen and e > A x 10~^ in the case 
of hydrogen. Therefore the reactions of Eqs. ([5^ are 
strictly in thermal equilibrium using the strictest of the 
two which is the one from hydrogen. In this case we can 
use the Saha equations that determines the equilibrium 
abundances 



ni 
n2 



A3 






'Eh/T ^ ^^ 



,-Eo/T ^ 5 



(33) 



n2 = 


- 25"^' 

"x 


ni 


1 + ^ + ^ 



where rii and 712 are the number densities of anoma- 
lous isotopes He and Oe respectively, and Aj:/(o) = 
^2-K jmiii^o'^T is the thermal De Broglie wavelength. 
From Eqs. (|33p one can solve for n\ and n^., 

A 



(34) 



where n^ is the total number density of millicharged par- 
ticles in the water (captured or free). There are strict 
constraints on ni/nn based on hydrogen anomalous iso- 
tope searches in sea water. These constraints set an up- 
per limit for ni/nn at 10^^^ — 10~^^ |5l|. However, 
it has been speculated that due to the heaviness of the 
milli-charged particles, anomalous isotopes might con- 
centrate at larger depths inside the oceans. Deep sea 
searches for hydrogen anomalous isotopes have imposed 
less strict upper limits of ~ 4 x 10~^^ [S^l, which are 
the ones we are going to use. If one uses Eqs. d^T]) . and 
(l34t . the ratio ni/nH can be estimated and compared 
to the upper limit from deep water searches of anoma- 
lous hydrogen isotopes. The result is quite interesting. 
For a TeV particle with relatively large e > 1.6 x lO"'' 
no constraints are imposed. Despite the relatively large 
values of e, due to the fact that the binding energy of 
the particle to the oxygen is much larger than that of 
the hydrogen (see Eq. (l27l) ). the millicharged particle 
strongly prefers to "hide" in oxygen rather than hydro- 
gen. Our estimate is in accordance with the qualitative 
argument presented in [17]. For the rest of the e phase 
space where the reactions (1321) are still in thermal equi- 
librium i.e. 4 X lO"'^ < e < 1-6 X lO""*, we find that TeV 
dark millicharged particles are excluded within the region 
2.9 X 10~s < e < 1-6 X lO^"*. Finally, below e < 4 x lO"'^, 
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there is no thermal equihbrium between the reactions 
(I32p . and therefore in that case the maximum value of 
ni/nu is given directly by Eq. (JSTJ (substituting p by 
Ph)- In reality photons or other type of collisions could 
destroy part of the He bonds reducing the ni abundance. 
However, even with no further distraction of the bound 
states no constraints can be imposed in this regime. The 
excluded region 2.9 x 10"^ < e < 1.6 x lO""* for a TeV 
particle, expands slightly to 1.8 x 10^^ ^ e ^ 1-7 x 10""^ 
for a 100 GeV millicharged particle. 

There are two caveats on the above derived exclusion 
region of e which might potentially lead to narrowing or 
even evasion of the constraints. The first one has to do 
with the nature of the experimental searches for anoma- 
lous isotopes. The examined water has gone through a 
phase of distillation and electrolysis among others, and 
therefore it is possible for the tiny bond of 10~^ eV (for 
an e ~ lO^'^, see Eq. ([271) ') between the millicharged par- 
ticle and the hydrogen to break apart. 

The second caveat has to do with the validity of 
Eq. ([29| . This is in principle valid as long as the Bohr 
radius of the system WIMP-nucleus is smaller than that 
of the electron-nucleus system. If the former Bohr ra- 
dius is larger than the latter, the capture cross section 
should be smaller, since it would be dominated by electric 
dipole type interactions between the millicharged parti- 
cle and the system nucleus-electron. The Bohr radius 
is a = aoTOe/(.^er7iAr)where niN is the mass of the nu- 
cleus and ao the usual atomic Bohr radius. It turns out 
that for hydrogen, an < oq if e > 5 x 10"* while for 
oxygen ao < ao if e > 4.2 x 10^^. If a > ao the cap- 
ture cross section is expected to be lower that that of 
Eq. (23) because among other, the binding energy of the 
millicharged-nucleus system will become smaller. The 
study of this type of dipole interactions goes beyond the 
scope of our paper. 

In our case of interest, searches for anomalous isotopes 
of oxygen might in principle give stricter constraints due 
to the fact that most of the millicharged particles will 
bind with oxygen in water. Although upper bounds on 
n^/no — 10^^^ exist in literature [53[, it is questionable 
if they can apply in our case due to the fact that the 
search took place within a limited range of masses and 
also due to the fact that the search did not take place 
on deep ocean waters. This is crucial since it is expected 
that i?20e will sink to large depths due to the heaviness 
of the millicharged particle. In addition, both caveats we 
presented for the hydrogen constraints apply also here. 
Even if one ignores the aforementioned issues about the 
deep water, mass range and breaking of the bond, we 
find for a millicharged mass of 100 GeV (1 TeV) that for 
e < 5 X 10^^ (e < 7.9 x 10^®) the oxygen constraints are 
evaded. 

Finally we examine possible constraints that can 
emerge from searches of anomalous isotopes of helium 
in the atmosphere. Searches for anomalous isotopes of 
helium in the air have also imposed strict upper limits 
on the ratio of concentrations of anomalous helium over 



the normal one thc of 10^^^ to 10~^^ 5J]. Following 
our treatment of the hydrogen and oxygen constraints 
from water, we can use Eq. (|24l) to estimate above what 
values of e millicharged particles slow down to thermal 
velocities. For 100 GeV (1 TeV) particle, this happens if 
e > 3.7 X 10"^ (e > 1.2 x 10"^). We examine the three 
basic reactions of interest 



He + e o Hee 
+ e^Oe 

N + e^ Ne. 



(35) 



since the two dominant components of atmospheric air is 
nitrogen, and oxygen. As before, one can estimate ithorm 
to establish thermal equilibrium among the reactions 
((35)) . For helium (which obviously gives the strictest time 
scale), in order for tthcrm << ^acc, e>5xl0~^. We have 
used an "effective height" for the atmosphere h = 20 km 
in order to get an estimate on uhc ■ Helium is significa ntly 
younger in the atmosphere with iacc = 2 x lO^years [54| . 
The corresponding Saha equations are 



ni 



^ie 



-EhJT ^ ^ 



A3 



-Eo/T ^ 



= S, 



A3 '^ - '-'' 



(36) 



where ni, n2, and n^ are respectively the densities of iJee, 
Oe, and Ne. The probability of capturing a millicharged 
particle into a bound state with an atmospheric nucleus 
inside the atmosphere is p = nah where n is the number 
density of the nucleus (oxygen, nitrogen or helium) in the 
atmosphere, a the cross section given by Eq. (I29|) and h 
is the effective atmosphere height. For example in the 
case of helium, riHe = Nne/^T^R^h, where Nhc is the 
total number of helium atoms in the atmosphere and R 
the radius of the earth. The total accumulated number 
of millicharged particles over Nhc is 

THe = {px/m)ATrR'^ti,ccVp/NHe = {px/m)t^ccV(T. (37) 

For e> 5 X 10^^ (i.e. if the reactions ((35t are in thermal 
equilibrium) , one can estimate the production of anoma- 
lous helium 



ni = 



1 + "/2A/B + 73A/C + A/uHe 



(38) 



where n^ is the total density of millicharged dark matter 
in the air. Despite oxygen representing a fraction 0.21 
(instead of 0.78 for nitrogen) of the atmospheric gas, it 
gives the highest contribution due to the Z dependence 
of Eq. dMl). The factors 72 = 2 x 0.21/(5.2 x lO"'^) 
and 73 = 2 X 0.78/(5.2 x 10^^) represent respectively 
the ratio of abundances of atomic oxygen over helium 
and atomic nitrogen over helium. The factors of 2 ac- 
count for the fact that oxygen and nitrogen are in the 



molecular form O2 and N2. Helium consists only a tiny 
5.2 X 10^® portion of the atmospheric air. The region 
6.8 X 10^5 < e < 1-5 X 10"'' where Eq. §Q are valid, 
is excluded. In this regime, the millicharged particles 
are captured efficiently by oxygen and nitrogen atoms 
and due to thermal equilibrium a significant amount is 
transferred to helium. As in the case of water, once e 
becomes sufficiently large (e > 1.5 x 10""'), due to the 
large binding energies of the millicharged particle with 
oxygen and nitrogen, no constraints from helium arise. 
For e < 5 x 10~^, although Eqs. (|36|) are not valid over- 
all, because reactions p5p are not in thermal equilib- 
rium, collisions of the Hee system with photons or other 
atoms might lead to destruction of the Hee population. 
In this case the first equation of ([36]) might give an up- 
per bound for ni. ni/nne can be easily estimated from 
Eq. (l38l) if one sets 72 = 73 = 0, and uses n^/nne = 'The 
from Eq. p7p . Upon making the above assumption no 
constraints can be imposed below e ^ 5 x 10~^. There- 
fore the exclusion region for a TeV millicharged parti- 
cle is 6.8 X 10"^ < e < 1-5 X 10""*. This expands to 
5 X 10-^ < e < 1-5 X 10-"* for a 100 GeV particle. 



V. CONCLUSIONS 

We presented a model of composite dark matter with 
millicharged constituents. We showed under what con- 
ditions such a candidate can provide a monochromatic 
130 GeV photon line via decays. We found that for a 
compositeness scale of order TeV, e ~ 10"^^. We also 
estimated that for e ^ 10~^, millicharged dark mat- 



ter can be produced thermally via the freeze-in mech- 
anism. Due to self-interactions this dark matter candi- 
date might be warmer than cold dark matter. We studied 
all possible constraints coming from self-interactions, the 
early Universe and BBN. We also studied possible con- 
straints from direct search experiments. Unlike what was 
thought until now, we demonstrated that millicharged 
dark matter cannot be excluded easily in underground 
detectors. For very small values of e the cross section 
is low, but even for intermediate values where the cross 
section can be sufficiently large, millicharged particles de- 
celerate substantially providing recoil energies below the 
threshold of the detectors. Finally we examined possible 
constraints arising from searches of anomalous hydrogen 
in deep ocean water, oxygen and helium in the atmo- 
spheric air. We found that searches of anomalous hydro- 
gen can potentially exclude the region (for a TeV particle) 
2.9 X 10~^ ^ e ^ 1-6 X 10"'*, while searches of anoma- 
lous helium could exclude a tiny region around e ~ 10~^. 
The candidate that provides the Fermi line, passes all 
these constraints. Millicharged dark matter thermally 
produced via the freeze-in mechanism seems to fall into 
the exclusion region of the hydrogen. However, these 
constraints have to be taken with caution since the ex- 
perimental procedure for the detection of the anomalous 
atoms might destroy the weak bond of the nucleus- WIMP 
system. In addition TeV millicharged dark matter with 
e ^ 10~^ can evade these constraints because it falls into 
the region which galactic magnetic fields prevent the par- 
ticles from entering the galactic disk and those initially 
trapped in the disc are accelerated and injected out of it 
relatively fast. 
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